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Abstract 

We study the quark energy loss in an expanding quark-gluon plasma. The expanding 
plasma produced in high energy AA-collision is described by Bjorken's model. The de- 
pendence of the energy loss on the infrared cutoff for the radiated gluons, on the quark 
mass, and on the initial conditions of QCD plasma is investigated. 



Study of the induced gluon radiation from a fast quark in a quark-gluon plasma (QGP) 
is of great importance in connection with the forthcoming experiments on high energy 
AA-collisions at the RHIC and LHC It is expected that the energy loss of high-p^ jets 
produced at the initial stage of ^4y4-collision may be an important potential probe for 
formation of QGP. 1 ' 2 -* In the recent works 3 ' 4 - 1 the quark energy loss, AE q , was estimated 
for a homogeneous QGP. Here we report on the evaluation of AE q in expanding QGP in 
the light-cone path integral approach to the induced radiation. 5 ) As in previous works, 
QGP is modelled by a system of static Debye screened scattering centers. ^ 

We consider a fast quark produced in the central rapidity region with a velocity perpen- 
dicular to the axis of yl^-collision. We choose the z axis along the initial quark velocity, 
and the quark production point is assumed to be at z = 0. Then, the distance passed by 
the quark in QGP, L = z, is close to the expansion time, r. For the r(2;)-dependence of 
the temperature of QGP we use prediction of Bjorken's model 6 -* Tr 1//3 = T t^ 3 . 

The probability of radiation of a gluon with the fractional longitudinal momentum x 
from a fast quark produced at z = is given by 5 - 1 
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Here K is the Green's function for the Hamiltonian (acting in the transverse plane) 

v( P ,z) = -i n< - z) " 3 [ P - X ' Z \ (3) 

and K v is the Green's function for the Hamiltonian (fj) with v(p, z) = 0. In (^|) the 
Schrodinger mass is fi(x) = E q x(l — x), Lf = 2E q x(l — x)/[m 2 q x 2 + m 2 g {l — x)] is the 
gluon formation length, here m q is the quark mass, and m g is the mass of the radiated 
gluon. We introduce the gluon mass to remove the contribution of the unphysical long- 
wave gluon excitations. In (g) n(z) is the number density of QGP, and <r 3 is the cross 
section of interaction of color singlet qqg system with color center. Summation over 
triplet (quark) and octet (gluon) color states is implicit in (H). The z-dependence of 03 is 
connected with the one of the Debye screening mass. 

The vertex factor g(z 1 ,z 2 ,x), entering ([I]), reads (we neglect the spin-flip q — > qg 
transitions, which give a small contribution to the quark energy loss) 
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The Hamiltonian ([!]) describes evolution of the light-cone wave function of a fictitious 
qqg system. 5 ) Using relations between the transverse Green's functions and the light-cone 
functions for the transition q — » qg in vacuum and inside medium 7 ) the radiation rate ([I]) 
after some algebra can be represented in another form 

lx- = J dZn{z) dx ' (5) 
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— ^j—^ — = Re J dpip*(p,x)a 3 (p,x,z)?p(p,x,z) , (6) 

where ip(p, x) is the light-cone wave function for the q — > qg transition in vacuum, 
and tp(p,x,z) is the medium-modified light-cone wave function for q —>■ qg transition 
in medium at the longitudinal coordinate z. In the low density limit and z — * oo 
ifi(p,x,z) becomes close to ip(p,x), and @ goes over into formula for the Bethe-Heitler 
cross section for a quark incident on an isolated center, da BH (x)/dx. 8 ^ At z — > the ra- 
tio dafM(x, z)/dx/da BH \x)/dx vanishes. 4 ' 8 ) This is a direct consequence of the decrease 
at small z of the transverse size of the qg Fock component of the quark produced at 
z = 0. Note that this effect is responsible for the L 2 - dependence of AE q at E q — > oo for 
a homogeneous medium. 3 ' 4 ) 

The three-body cross section entering the imaginary potential (|3|) can be written as 9 -* 

9 1 

a 3 (p,x,z) = ~[a 2 (p,z) + «7 2 ((1 - x)p, z)\ - -a 2 (xp,z), (7) 

where a 2 (p,z) is the dipole cross section of interaction with color center of color singlet 
qq system. The latter can be written as a 2 (p,z) = C 2 (p,z)p 2 , where C 2 (p,z) has a 
smooth (logarithmic) dependence on p at small p. As in Ref. 4 we approximate C 2 (p,z) 
by its value at p ~ l/m g . Then the Hamiltonian (^j takes the oscillator form with the 
^-dependent frequency 
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where C 3 (x, z) — | {9[1 + (1 — x) 2 } — x 2 } C 2 (l/m g , z) . The coefficient C 2 was calculated 
in the double gluon approximation. In numerical calculations we evaluated the Green's 
function K using the approach previously developed in analysis of the vector meson 
photoproduction. 10 ) 

It follows from (pj) that even for a thermolized QGP the oscillator frequency is dom- 
inated by the gluon contribution. There is every indication 11,12 ) that for the RICH and 
LHC the hot QCD medium produced in AA-collisions at r ~ r ~ 0.1 fm will be ther- 
molized gluon plasma to a first approximation, and for quarks the chemical equilibration 
is not reached during the expansion of QGP. For this reason we take for the gluon fugacity 
X g = 1. To take into account the suppression of quarks we use the value X q = 1/3 for 
the quark fugacity. Numerical calculations were carried out with a s = 1/3. For tq we 
use the value 0.1 fm. 11,12 ) At z < r we take Q(z) = f2(r ). We have evaluated AE q for 
Tq = 1100 expected for central Pb-Pb collisions at the LHC. 11 - 1 To study the dependence 
of the quark energy loss on To we give also the results for To = 700 MeV. The numerical 
predictions for AE q as a function of E q for L = 3, 6, 9 fm are shown in Fig. 1. The results 
were obtained for m q = 0.2 GeV. As for a homogeneous QGP 4 \ our predictions have a 
weak dependence on m q . To study the infrared sensitivity of AE q we present in Fig. 1 the 
results for m g = 0.75 and 0.375 GeV. These values are of the order of the Debye screening 
mass for QGP in the region r > 2 — 3 fm, which, as our numerical calculations show, 
dominates the quark energy loss. For this reason the value of m g within the above range 
seems to be a plausible estimate for the infrared cutoff in the considered problem. Note 



3 



that in the limit E q — > oo AE q has only a smooth (logarithmic) m 9 -dependence connected 
with p-dependence of the coefficient C^. 4 - 1 Our numerical calculations really demonstrate 
that for E q > 10 GeV AE q is not very sensitive to m g . Fig. 1 shows that AE q grows by 
a factor ~ 2 — 3 as E q increases from ~ 10 to ~ 100 GeV. The predictions for L = 6 and 
L = 9 fm can be regarded as a plausible estimates for AE q in central collisions of heavy 
nuclei. We have checked that the effect of the mixed phase (with T c « 150 MeV) and the 
hadronic phase turns out to be relatively small. Comparison of the results for T = 1100 
and T = 700 MeV demonstrates that the T -dependence of AE q is not very strong. This 
is connected with an increase of Landau- Pomeranchuk-Migdal suppression and a decrease 
of the coefficient C 2 for a high density QGP. 

To study the sensitivity of AE q to dynamics of QGP at small r we have also carried 
out the calculations taking Q(z < a) = Q(a) with a ~ 1 — 2 fm. The results for these 
versions do not differ strongly from the ones in Fig. 1. This demonstrates that AE q 
is insensitive to the dynamics of QGP at r < 2 fm. This fact is a consequence of the 
vanishing of the ratio da^{x, z) / dx / da BH \x)/dx at small z. 

To illustrate the quark mass dependence of AE q in Fig. 2 we compare the results for 
AE q for c-quark (m q = 1.5 GeV) with the ones for light quarks (m q = 0.2 GeV). As one 
can see the dependence of AE q on m q becomes weak at E q > 50 — 100 GeV. 

In Fig. 3 we present the gluon spectra for m q = 0.2 and 1.5 GeV at E q =50, 100 and 
200 GeV for T = 1100 MeV. We also show in this figure the Bethe-Heitler spectra. As 
one can see the radiation rate is strongly suppressed by the medium and finite-size effects 
for the gluon momenta k > 10 GeV (excluding a narrow region near k ~ E q ). For the 
heavy quark the radiation rate at k « E q is suppressed as compared to the light quark. 
However, in the soft region the spectra for the light and heavy quarks are close to each 
other. 

Our predictions for the absolute normalization of the energy loss must be regarded 
as rough estimates with uncertainties of a factor ~ 2. Nevertheless, the obtained rather 
large values of AE q show that jet quenching may be an important probe for the formation 
of QGP in AA-collisions. The Monte-Carlo analysis of jet quenching 2 ) with AE q close to 
the estimates of the present work demonstrates that the induced radiation must consid- 
erably modify the charged particle spectra. The induced gluon radiation can also lead 
to other interesting experimental consequences. For instance, fluctuations of the quark 
energy loss will generate additional transverse momentum (defined with respect to the 
AA-collision axis) in production of cc and bb pairs. In the case of g — > gg transitions 
the production of the gg state in the color decuplet state can increase the cross section 
for production of baryon-antibaryon pairs through the mechanism analogous to the one 
previously discussed 13 - 1 in connection with BB annihilation at high energies. 

This work was partially supported by the INTAS grants 93-239ext and 96-0597. 
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Figures 
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Figure 1: The quark energy loss as a function of the quark energy for m q = 0.2 GeV, 
m g = 0.75 GeV (solid line) and m g = 0.375 GeV (dashed line). 
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Figure 2: The quark energy loss as a function of the quark energy for m g = 0.75 GeV, 
m q — 1.5 GeV (solid line) and m q = 0.2 GeV (dashed line). 
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Figure 3: The gluon spectrum (solid line) as a function of the gluon momentum at E q =50, 
100 and 200 GeV for m q = 0.2 (a) and 1.5 (b) GeV. In both the cases m g = 0.75 GeV, 
and Tq = 1100 MeV. The dashed line shows the Bethe-Heitler spectrum. 
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